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CONVERSION  FACTORS,  U.  S.  CUSTOMARY  TO  METRIC  (SI) 
UNITS  OF  MEASUREMENT 


U.  S.  customary  units  of  measurement  used  in  this  report  can  be  converted 
to  metric  (SI)  units  as  follows: 


Multiply 

By 

To  Obtain 

Acres 

4046.873 

Square  meters 

Cubic  feet  per  day  per 
foot 

0.0929 

Cubic  meters  per  day  per 
meter 

Feet 

0.3048 

Meters 

Feet  for  second 

0.0348 

Meters  per  second 

Horsepower  (550  foot¬ 
pounds  per  second) 

745.6999 

Watts 

Miles  per  hour 
(U.  S.  statute) 

1.609347 

Kilometers  per  hour 

Square  feet 

0.09290304 

Square  meters 

Square  feet  per  second 

0.09290304 

Square  meters  per  second 

COMPUTER  MODELING  OF  HYDRODYNAMICS  AND  SOLUTE  TRANSPORT 
IN  CANALS  AND  MARINAS 


Literature  Review  and  Guidelines  for  Future  Development 
I.  INTRODUCTION 


1. 1  Background 

Over  the  past  few  decades,  there  has  been  a  tremendous  increase  in  the  de¬ 
mand  for  residential  canals  and  marinas  along  the  coasts  of  the  United 
States.  This  demand  has  led  to  the  design  and  construction  of  numerous  such 
systems,  often  without  the  benefit  of  a  sound  scientific  background. 

As  an  example,  from  the  late  1940's  to  the  early  1960 1  s ,  many  large-scale 
residential  waterfront  developments  were  constructed  in  Florida.  The  great 
majority  of  these  developments  were  branching  networks  of  canals  (Figure 
1.1).  They  were  built  by  dredging  and  filling  coastal  wetlands.  By  1974, 
two  developments  of  over  100,000  acres*  each,  and  14  over  10,000  acres, 
were  either  occupied  or  under  construction.  The  total  land  area  under  devel¬ 
opment  in  Florida  alone,  in  1974,  exceeded  810,000  acres  (Carter,  1974). 

Unfortunately,  it  is  now  evident  that  many  of  these  systems  are  not  the 
beautiful  water  bodies  anticipated,  but  rather  they  are  beset  with  all 
kinds  of  problems.  Many  of  these  systems  flush  poorly,  have  poor  water 
quality,  are  silting  up,  or  have  weak  circulation.  The  same  people  who 
demanded  that  such  systems  be  built  are  now  demanding  that  scientists  and 
engineers  study  them  in  more  detail  to  learn  how  they  work  and  what  consti¬ 
tutes  a  good  design. 


*  A  table  of  factors  for  converting  U.S.  customary  units  of  measurement  to 
metric  (SI)  is  presented  on  page  vi. 


1-1 


1-2 


Up  to  this  time,  the  "science"  of  canal  and  marina  design  in  the  coastal 
zone  has  grown  largely  by  learning  from  past  mistakes.  Analytic  and  numer¬ 
ical  models  used  were  often  developed  for  other  purposes,  and  were  not 
always  adequate  to  describe  the  physics.  The  contractor's  prime  motivation 
was  to  maximize  the  number  of  waterfront  dwellings  or  boat  slips. 

The  problem  is  further  compounded  by  the  permitting  process,  in  which 
local,  State,  and  Federal  agencies  must  evaluate  canal  and  marina  designs 
and  the  validity  of  analytic  tools  used.  Given  a  lack  of  understanding  of 
the  physics  of  these  systems,  and  the  multitude  of  "simulation"  tools  that 
can  and  have  been  used,  it  is  no  wonder  that  the  permitting  agencies  are 
somewhat  confused,  and  in  some  cases  resort  to  a  position  of  "no  more"  as  a 
way  out.  This  is  a  safe  approach  because  it  is  easy  to  suggest  that  we  will 
negatively  impact  the  environment,  as  we  have  in  the  past. 

The  design  engineer's  view,  however,  is  that  they  can  now  build  better  sys¬ 
tems  that  are  compatible  with  the  environment,  and  indeed  in  some  cases  may 
improve  on  it. 

Consequently,  within  the  past  few  years,  much  more  emphasis  has  been  placed 
on  learning  what  makes  canal  and  marina  systems  work.  Although  the  surface 
has  only  just  been  scratched,  today  there  is  a  better  understanding  of  the 
physics  involved,  and  thus  a  better  understanding  of  how  to  develop  the 
analytical  tools  needed  for  canal  and  marina  design. 

I .2  Purpose  and  Scope 

The  Corps  of  Engineers,  as  a  permitting  agency,  must  be  in  an  informed 
position  to  evaluate  proposed  canal  and  marina  designs.  An  understanding  of 
the  physical  processes  at  work,  and  of  the  ability  of  existing  analytical 
tools  to  accurately  simulate  these  processes,  is  essential.  Such  judgments 
must  be  made  in  the  light  of  local,  State,  and  Federal  regulations.  The 
discussion  of  present  conditions  is  covered  in  Chapters  1 1 -I  V . 
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One  of  the  purposes  of  this  report  is  to  discuss  the  present  understanding 
of  how  these  systems  work,  and  then  review  existing  numerical  and  analyti¬ 
cal  techniques  for  simulating  hydrodynamic  and  solute  transport  processes. 
To  do  this,  past  studies  of  canal  and  marina  designs  and  simulations  will 
be  reviewed,  and  such  evaluation  criteria  as  flushing  rates  and  residence 
time  will  be  examined. 

As  a  research  and  modeling  group,  the  U.S.  Army  Engineer  Waterways  Experi¬ 
ment  Station  (WES)  is  in  a  position  to  extend  our  ability  to  understand  and 
model  these  systems.  In  this  report,  future  needs  in  canal  and  marina  de¬ 
sign  and  simulation  will  be  addressed.  This,  with  an  approach  to  further 
studies,  will  be  the  subject  of  Chapter  V. 

I .3  Distinction  Between  Canals  and  Marinas 

In  this  report,  a  distinction  will  be  made  between  the  definitions  of 
canals  and  marinas,  although  it  is  realized  that  such  a  difference  may  not 
always  be  significant.  This  is  done  based  on  geometric  and  usage  consider¬ 
ations. 

Christensen  and  Snyder  (1978)  have  classified  the  above  types  of  water  bod¬ 
ies  into  seven  major  groups  (Figure  1.2): 

1.  Simple  dead-end  canal, 

2.  Flow-through  canal, 

3.  Comb- structure  canal  system, 

4.  Higher  order  finger  canal, 

5.  Canal  with  lagoon  basin, 

6.  Lagoon  with  two  tidal  entrances,  and 

7.  Borrow  lake. 

A  canal  is  considered  to  be  a  long,  narrow  waterway,  which  is  primarily 
used  for  residence  with  a  docking  facility.  Such  canal  systems  (numbers  1-4 
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DEAD-END  CANAL 


INTRACOASTAL  WATERWAY 
0* 

OCEAN 


FLOlWHROUGH  CANAL 


Figure  1.2.  Classification  of  Water  Bodies 
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above)  are  frequently  found  along  the  lower  Atlantic  and  Gulf  Coast  sea¬ 
boards  where  the  tidal  ranges  are  small.  Waterfront  canals  will  be  further 
subdivided  into  three  types  (Barada  and  Partington,  1972;  Lindall  and 
Trent,  1975): 

1.  Bay-fill  or  finger  fill  -  those  constructed  below  mean  low  tide  by 
dredging  and  filling  shallow  bay  bottoms  (Figure  I. 3a), 

2.  Intertidal  -  those  constructed  by  dredge-and-fill  between  mean  low  and 
mean  high  waters.  In  many  cases,  these  canals  are  located  in  mangrove 
or  salt  marsh  ecosystems,  in  bays,  estuaries,  lakes,  or  other  wetlands 
(Figure  I .3b) ,  and 

3.  Inland  or  upland  -  those  developed  by  excavating,  which  are  above  mean 
high  water  and  connected  to  natural  channels,  lakes,  rivers,  or  other 
natural  or  artificial  waterways  (Figure  I. 3c). 

Marinas,  on  the  other  hand,  are  primarily  used  for  boat  docking.  The  width 
is  often  on  the  same  order  of  magnitude  as  the  length  (Figure  1.5).  Classi¬ 
fications  5  and  6  in  Figure  1.2  show  typical  shapes. 
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c.  Inland  Development 
Figure  1.3.  Waterfront  Canals 
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Figure  1.4.  Schematic  Examples  of  Pacific  Northwest  Marina  Geometries 
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II.  PROCESSES 


1 1 .1  Canals 

The  flows  in  residential  canal  systems  are  complex,  time-varying  interac¬ 
tions  of  forcing  due  to  tide,  wind,  density-induced  currents,  and  secondary 
currents  due  to  topographic  features  such  as  bends  and  sills.  An  under¬ 
standing  of  these  processes,  and  their  relative  effects,  is  crucial  in 
choosing  both  design  approaches  and  evaluation  criteria. 

Several  different  institutions  have  made  fairly  comprehensive  measurements 
in  canal  systems,  although  such  studies  tend  to  concentrate  on  one  or  more 
aspects  of  canal  processes  and  omit  others.  Considering  hydrodynamics  and 
transport  mechanisms,  a  number  of  measurements  were  performed  by  the  Hy¬ 
draulic  Laboratory  of  the  University  of  Florida  during  a  three-year  Sea 
Grant  and  other  hydrographic  studies  (University  of  Florida,  1976;  Morris 
and  Christensen,  1976;  Morris  et  al.,  1977,  1978a;  Walton  et  a 1 . ,  1975a, 
1975b).  Typical  canal  parameters  are  listed  in  Table  II. 1. 

Using  tidal  range  as  a  measure  of  available  energy  to  these  systems,  resi¬ 
dential  canals  along  the  Florida  and  Gulf  Coasts  are  classified  as  low 
energy  canals.  Tidal  ranges  (Figure  II. 1)  are  small,  as  are  induced  tidal 
velocities,  even  at  the  mouths  of  many  systems.  Maximum  velocities  on  the 
order  of  0.5  ft/sec  are  not  uncommon,  and  these  magnitudes  decrease  towards 
the  dead  ends  of  these  systems.  Observed  maximum  water  surface  slopes  were 

_5 

measured  to  be  on  the  order  of  10  . 

Because  induced  tidal  velocities  are  relatively  small,  other  types  of  forc¬ 
ing  can  often  play  a  significant  role  in  canal  hydrodynamics.  In  many  of 
the  systems  observed  (Walton  et  al.,  1975b,  Morris  et  al.,  1978a),  wind 
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Table  II. 1  -  Typical  Measured  Canal  Parameters 


PARAMETER 

RANGE 

Length,  L 

1,500  -  10,000  ft 

Bottom  width,  b 

50  -  100  ft 

Mean  tidal  depth,  dQ 

5  -  12  ft 

Inverse  side  slope,  s 

0-5 

Tidal  range,  r:  Atlantic 

2  -  5  ft 

Gulf  of  Mexico 

2  -  3  ft 

Tidal  period,  T 

12.42  hr 

Nikuradse's  equivalent  sand  roughness,  k 

1  -  20  ft 

Maximum  (in  canal )  longitudinal  dispersion 
coefficient,  E 

0.5  -  5.0  ft2/sec 

Dimensionless  dispersion  coefficient,  K 

2  -  20 

Maximum  water  surface  slope,  S 

10"5  -  10~6 

Table  II. 2  -  Typical 

Measured  Marina  Parameters 

PARAMETER 

RANGE 

Length,  L 

1000  - 

3000  ft 

Width,  b 

400  - 

1000  ft 

Mean  lower  low  water  depth 

5  - 

15  ft 

Tidal  range,  r 

3  - 

13  ft 

Tidal  period,  T 


12.4  hr 


ELEVATION  (feet) 


>1*.  Fla  Octabar  74  -  77.  1975 


Muni,  Fla  Ociaba*  74-77.  1*75 


S'.  FinriSa'f,  Fla  Ociabar  74  -77.  1975 


TIME  (hours) 

Figure  II. 1.  Typical  Tide  Curves 
in  Floridian  Canals 
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was  perhaps  the  dominant  short-term  driver.  When  a  wind  blows  along  the 
axis  of  a  dead-end  canal  system  for  some  period  of  time,  continuity  of  mass 
may  produce  a  counterflow  in  the  lower  layers  (Figure  11.2).  In  a  flow¬ 
through  system,  the  wind  will  tend  to  superimpose  on  the  tidal  flow,  in  the 
one-dimensional  sense  of  longitudinal  transport. 

This  observation  of  wind-induced  counterflows  may  be  seen  in  other  measure¬ 
ments  as  well.  Dye  introduced  into  a  canal  system  on  a  windy  day  may  not 
exhibit  the  Gaussian  distribution  expected  of  one-dimensional,  mean  flow 
physics  (Figure  1 1.3).  In  this  case,  the  wind  was  blowing  down  the  canal 
producing  an  observable  flow  reversal  in  the  lower  layers  (Figure  11.2). 
This  carried  some  of  the  dye  upstream  where  it  was  upwelled  or  diffused 
into  the  surface  water  to  elongate  the  "tail"  of  the  distribution.  This 
same  physical  explanation  can  be  given  to  explain  the  data  of  dye  studies 
performed  by  the  Environmental  Protection  Agency  (EPA  1975). 

A  longer  term  mechanism,  observed  particularly  in  south  Florida  canals  dur¬ 
ing  the  hot,  wet  summer  months,  is  density- induced  circulation.  McKeehan 
(1975),  van  de  Kreeke  et  al .  (1977),  and  Monopolis  (1978)  have  noticed  that 
during  these  periods  of  relatively  light  winds  and  relatively  large  fresh¬ 
water  inflows,  density  stratification  can  occur  (Figure  II. 4),  which 
induces  circulation.  Van  de  Kreeke*  believes  that  these  density- induced 
circulations  are  important  in  long-term  circulation  and  water  quality 
transport,  particularly  during  the  wet  surnner  months. 

Even  superimposed,  these  various  types  of  forcing  produce  real  time  veloci¬ 
ties  whose  magnitudes  are  generally  small.  Examples  of  three-layer  flows  in 
these  shallow  systems  are  not  uncommon  (Figure  II. 4).  These  small  veloci¬ 
ties  in  turn  lead  to  relatively  low  rates  of  mixing  as  seen  from  observed 
longitudinal  dispersion  coefficients.  Dye  studies  in  a  number  of  canals 

(Morris  et  al . ,  1978a;  Walton  et  al . ,  1975a,  1975b;  EPA,  1973,  1975)  have 

2 

measured  coefficients  on  the  order  of  5  ft  /sec  and  less  in  these  systems. 


*  Personal  Communication,  J.  van  de  Kreeke,  Rosenstiel  School  of  Marine 
and  Atmospheric  Science,  University  of  Miami,  Coral  Gales,  FL,  1982. 
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Figure  II. 2.  Salinity  and  Velocity  Measurements  Show  Flow 
Reversal  in  Loxahatchee  Canal,  Florida 


Figure  II. 3.  Non-Caussian  Distribution  of  Dye 
Concentration  in  Floridian  Canal 
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This  relatively  weak  mixing  due  to  variations  in  the  cross-sectioned  velo¬ 
city  distribution  means  that  understanding  the  induced  mixing  effects  of 
topographic  features  such  as  bends,  boat  docks,  branches,  and  undulating 
beds  becomes  important.  A  number  of  laboratory  studies  at  the  University 
of  Florida  (O'Hargan,  1975;  Parr  and  Christensen,  1977;  Swenty,  1977)  have 
shown  how  longitudinal  dispersion  may  be  increased  by  the  presence  of  bends 
and  boat  docks.  This  effect  has  also  been  noted  in  streams  (Fischer,  1967, 
1969). 

It  is  interesting  to  note  on  this  subject  that,  even  though  these  systems 
are  classified  as  low  energy  tidal  canals,  the  available  energy  used  in 
friction  and  other  losses  is  small,  as  can  be  calculated  using  classical 
hydraulic  formulae.  Defining  the  potential  energy,  E  ,  of  the  tides  in  a 
canal  system: 


Ep  =  As  y  a2  (II. 1) 

2 

where  A  =  area  of  water  surface  (L  ), 

s  2 

y  =  weight  density  of  water  (ML/T  ),  and 

a  =  tidal  amplitude  (L). 

Snyder*  calculated  the  available  energy  to  a  canal  system  he  was  studying 
to  be  1.4  hp.  He  further  calculated  that  losses  due  to  form  drag,  bends, 
constrictions/expansions,  pipes/culverts,  and  friction  were  only  0.5  per¬ 
cent  of  the  available  energy. 

In  a  series  of  experiments,  Morris  et  al .  (1978a)  deployed  three  dual-axis 
electromagnetic  current  meters  in  a  cross  section  of  a  canal  just  below  a 
sharp  bend.  Periodically,  vertical  profiles  were  taken  at  each  station  and 
pieced  together  to  produce  a  cross-sectional  view  of  the  velocity  field 
(Figure  1 1.5).  From  this  figure,  we  can  clearly  see  the  induced  secondary 
current,  which,  if  persistent,  would  overturn  the  water  column,  bringing 
oxygen-depleted  waters  to  the  surface. 


*  Memorandum  to  Florida  Department  of  Environmental  Regulation  from  Snyder 
Oceanography  Services,  Jupiter,  FL,  17  June  1978. 
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Figure  II. 5.  Secondary  Current  in  57  Acres  Canal 
System,  Florida 
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A  knowledge  of  canal  hydrodynamics  is  important  not  only  for  itself,  but 
also  in  understanding  water  quality  and  sedimentation  processes.  The  liter¬ 
ature  on  the  environmental  effect  of  canals  and  dredged  holes  contains  data 
on  salinity,  temperature,  and  dissolved  oxygen  which  are  useful  for  the  de¬ 
sign  engineer.  Some  of  the  literature  was  first  reviewed  for  a  study  of  the 
environmental  impact  of  borrow  pits  in  Maryland  estuarine  waters  (Poll's, 
1974).  Pol  is  summarized  work  in  Texas  and  on  the  Atlantic  Coast  from  Flor¬ 
ida  to  New  Jersey.  Since  then,  much  additional  work  has  been  performed  on 
studying  water  quality  in  tidal  canal  systems.  The  EPA  (1973,  1975)  did 
perhaps  the  largest  study  in  canals  from  Florida  to  North  Carolina.  Paulson 
et  al .  (1974,  1975,  1977a,  1977b)  measured  water  quality  parameters  in 
canals  in  Mississippi,  and  Bailey  (1977)  did  some  recent  work  in  Florida. 

Dissolved  oxgyen  (DO)  is  often  considered  to  be  a  standard  of  canal  health, 
and  values  anywhere  between  zero  and  saturation  have  been  measured.  Zero 
values  are  often  associated  with  the  dead  ends  of  canal  systems  and  deep 
holes  (Lindall  et  al . ,  1975  -  Figure  II. 6).  Low  DO  values  have  been  held 
responsible  for  fish  kills  in  Floridian  canals  where  it  is  thought  that  the 
fish  become  trapped  and  cannot  find  their  way  to  oxygen  rich  waters. 

The  most  comprehensive  set  of  DO  measurements  in  canals  attempted  to  relate 
DO  values  to  the  depth  of  the  canal  system  to  determine  maximum  desired 
depths  in  these  systems  (EPA,  1973,  1975).  Results  of  this  analysis  are 
shown  in  Figures  II. 7  and  II. 8  for  North  Carolinian  and  Floridian  canals, 
but  as  can  be  seen  in  Figures  1 1. 9  and  1 1. 10,  for  two  different  canal  sys¬ 
tems  in  Florida,  this  result  is  not  universal.  Such  an  analysis  provides  a 
reasonable  basis  for  a  design  regulation.  However,  more  work  needs  to  be 
done  to  understand  why  some  systems  differ  from  the  general  observation, 
and  under  what  conditions. 

In  an  interesting  study  of  water  quality.  Messenger  and  Reynolds  (1979) 
obtained  five  months  of  data  in  a  residential  canal  network  on  the  Texas 
coast  (Figure  1 1 . 1 1 ) .  They  attempted  to  relate  measured  values  to  the  ages 
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Figure  II. 6.  Monthly  DO 
Readings  (from  Linda  11 
et  al.,  1975) 


Figure  II. 7.  DO  Values  in  North 
Carolina  Canals  Varying  With 
Depth  (from  EPA,  1975) 
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Figure  II. 8.  DO  Values  in  Florida  Canals  Varying 
With  Depth  (from  EPA,  1975) 
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Figure  II. 10.  Examples  of  DO  Profiles  Above  and  Below  Standards 
(4  mg/1)  in  Punta  Gorda  Canals,  Florida  (from  EPA,  1975) 


of  the  canals  studied  in  the  system.  Their  results  (Figure  11.12)  indicated 
that  age  was  not  a  dominant  factor  in  determining  canal  water  quality  when 
flushing  and  circulation  are  also  exerting  an  effect.  They  believed  this  to 
be  the  case  in  this  system,  where  the  older  canals  are  built  closer  to  the 
receiving  water  body,  and  generally  have  better  flushing  properties. 

Two  physical  features  that  have  received  much  attention  are  canal  depth  and 
the  presence  of  a  sill.  When  a  canal  is  first  dredged,  before  connection  to 
a  receiving  water  body,  a  plug  is  often  left  in  place.  Upon  removal ,  a  sill 
may  remain  which  impedes  the  circulation  of  the  bottom  waters.  As  a  result, 
sediment  may  be  trapped. 

It  has  been  observed  that  "deep"  canals  are  not  adequately  flushed  by  tidal 
action,  and  that  lower  layers  act  as  a  trap  for  sediments  and  organic  de¬ 
tritus.  Pol  is  (1974)  and  Barada  and  Partington  (1972)  reported  thermal 
stratification  in  canals  deeper  than  15  feet.  A  sharp  interface  was  mea¬ 
sured  between  10-12  feet,  with  indications  of  less  turbidity,  anaerobic 
conditions,  and  the  presence  of  hydrogen  sulfide  below  the  interface.  It 
has  also  been  observed  that  canals  that  are  very  shallow  (under  5  feet)  may 
have  poor  flushing  characteristics,  poor  navigability,  and  increased  tur¬ 
bidity  due  to  boat  traffic  (Chesher,  1974). 

Sedimentation  has  been  observed  in  many  tidal  canals.  Such  buildup  may  be 
caused  by  direct  stormwater  discharge  into  these  systems,  or  high  sediment 
loads  from  the  receiving  waters.  Sediment  buildup  is  usually  associated 
with  topographic  features  such  as  deep  holes,  sills,  and  bends,  all  of 
which  tend  to  slow  the  water.  On  one  occasion,  a  navigation  project  on  the 
Intracoastal  Waterway  in  Florida  caused  a  side  canal  to  silt  up  at  its 
mouth  and  be  impassible  at  low  tide  (Walton  et  al.,  1975a). 

The  hydrodynamics  and  mass  transport  properties  of  tidal  canal  networks  are 
closely  interrelated.  A  thorough  knowledge  of  these  processes  is  needed  to 
understand  how  to  build  these  systems  to  reduce  water  quality  problems  and 
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to  improve  flushing  and  turnover  where  needed.  A  knowledge  of  these 
processes  is  needed  if  we  are  to  be  able  to  develop  good  predictive  tools 
for  canal  design  and  evaluation. 


II .2  Marinas 

Many  of  the  physical  processes  observed  in  canal  systems  are  important  in 
marina  hydrodynamics  and  mass  transport.  Generally,  however,  time  scales 
may  differ  due  to  a  combination  of  geometry  and  tidal  ranges,  and  processes 
may  have  a  different  order  of  importance.  Typical  geometric  parameters  are 
shown  in  Table  II. 2. 

Marinas  are  usually  built  adjacent  to  large  water  bodies,  such  as  estuaries 
and  oceans,  to  provide  protection  and  calm  water  for  boat  launching  and 
storage.  The  circulation  between  the  marina  and  the  adjacent  water  body  is 
a  combination  of  tidal,  density,  and  current  effects.  These  processes  are 
described  by  Vollmers  (1976),  Westrich  (1976),  Askren  (1979), and  Nece  et 
al .  (1976)  and  are  shown  in  Figures  11.13  and  11.14. 

Tidal  ranges  associated  with  marina  locations  are  much  more  varied  than  for 

most  canal  systems.  They  vary  from  the  low  energy  systems  found  along  the 

lower  Atlantic  and  Gulf  Coasts  to  the  relatively  larger  ranges  found  along 
the  Pacific  Coast  and  in  New  England.  Ranges  in  the  Puget  Sound  marinas 
study  by  Nece  et  al .  (1972,  1974,  1975),  for  example,  were  between  3  and  13 
feet  depending  on  location.  These  relatively  larger  tidal  ranges,  together 
with  the  same  order  of  low  water  depths  for  navigation,  mean  that  flushing 
times  are  usually  smaller  (Figure  11.15). 

The  current  effect  (Figure  11.13)  is  not  often  noted  in  canal  systems  be¬ 
cause  they  are  relatively  narrow.  However,  in  many  marinas  there  may  exist 

a  velocity  shear  between  the  basin  and  the  adjacent  water  body  which  may 

produce  one  or  more  circulation  cells  (vortices  -  Figure  11.16).  The  theo¬ 
retical  formation  of  such  a  circulation  cell  is  described  by  Askren  (1979) 
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Figure  11.13.  Components  of  Marina  Circulation 
(after  Vol lmers ,  1976) 


Figure  11.14.  Induced  Circulation  Pattern  in  Side  Basin 
(from  Askren,  1979;  and  Westrich,  1976) 


11-17 


TIDAL  RANK.  Fee T 


gure  11.15.  Flushing  Characteristics  of  Various  Puget 
Sound  Marinas  (from  Nece  et  al . ,  1975) 
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Figure  11.16.  Drogue  Pathlines  in  Edmonds  Marina 
(from  Nece  and  Knoll,  1974) 


and  Westrich  (1976),  and  is  shown  in  Figure  11.14.  Such  cells  can  have  an 
effect  on  water  quality  and  sedimentation  within  the  marina,  but  at  this 
point  their  physics  are  not  well  understood. 

Wind  effects  over  marinas  tend  to  be  relatively  unimportant,  except  in 
driving  the  adjacent  water  body.  Over  the  marina  itself,  fetch  lengths  are 
small,  and  there  is  generally  much  masking  due  to  trees  and  houses. 
Density- induced  circulations  are  also  small,  as  little  fresh  water  enters 
these  systems.  Marina  water  may  tend  to  be  slightly  warmer  than  that  of  the 
adjacent  water  body  (Askren,  1979),  but  generally  not  sufficiently  so  to 
generate  significant  currents  and  changes  in  flushing  rates  in  systems 
where  tidal  exchange  dominates  (Nece  et  al . ,  1979). 

As  marina  development  has  grown,  researchers  are  noticing  more  and  more 
that  the  design  which  allows  for  calm,  protected  waters  also  acts  as  a  sed¬ 
iment  and  pollutant  trap.  Ecologic  and  water  quality  studies  are  relatively 
few  (Bowerman  and  Chen,  1971;  Nece  and  Knoll,  1974;  Nece  et  al . ,  1975),  but 
results  indicate  that  longer  detention  times  can  impair  water  quality  (Fig¬ 
ure  11.17).  In  studies  in  West  Coast  marinas,  Slotta  and  Noble  (1977), 
Schluchter  and  Slotta  (1978),  and  Swartz*  have  noted  large  amounts  of  ben¬ 
thic  deposition,  perhaps  up  to  2  meters.  There  is  also  some  evidence  that 
toxins  may  be  building  up  in  these  systems,  but  this  is  a  subject  of  on¬ 
going  research.**  Ironically,  it  appears  that  benthic  deposition  may  not 
improve  with  "good"  flushing,  but  that  it  may  actually  increase  as  more 
sediment- ladened  waters  are  interchanged  with  the  marina  basin.  Further 
studies  have  also  shown  that  boat  props  may  tend  to  resuspend  benthic  mate¬ 
rial,  particularly  in  shallow  marinas,  and  move  materials  around  the  basin. 
Such  findings  are  preliminary,  and  are  the  first  interpretations  of  long¬ 
term  ongoing  research  into  these  and  related  processes. 


*  Personal  Communication,  R.  Swartz,  EPA,  Corallis,  OR,  1982. 

**  Personal  Communication,  L.S.  Slotta,  Oregon  State  University,  Corvallis, 
OR,  1982. 
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Figure  11.17.  Observations  of  Water  Quality  Parameters 
in  Lagoon  Point  Marina  (Nece  et  al.,  1975) 
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An  engineering  approach  to  overcoming  sediment  buildup  in  some  marinas  on 
the  East  Coast  is  to  design  the  basin  with  a  lock  structure.  The  object  is 
to  reduce  the  exchange  of  sediment- laden  water,  and  also  to  reduce  tidal 
fluctuations.*  If  this  is  done  in  combination  with  good  engineering  prac¬ 
tice  (van  de  Kreeke  and  Larsen,  1981),  such  as  limiting  pollutant  loadings 
to  the  basins  like  stormwater  runoff  and  allowing  flushing  using  an  up¬ 
stream  water  supply,  good  water  quality  and  low  sedimentation  rates  can 
result. 


*  Personal  Communication,  B.  Edge,  Cubit  Engineering,  Inc.,  Clemson,  SC, 
1982. 
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III.  DESIGN  CONSIDERATIONS 


Canal  and  marina  design  and  construction  must  be  accomplished  within  the 
limits  of  local.  State,  and  Federal  regulations.  It  is  not  our  purpose  to 
review  these  here,  as  regulations  vary  widely  between  states  and  even  with¬ 
in  states,  but  instead  to  discuss  design  considerations  in  the  context  of 
analytical  evaluation  of  proposed  systems. 

Design  standards  and  environmental  regulations  have  tended  to  develop 
through  a  process  of  "learning  by  our  mistakes."  In  some  cases,  reaction  to 
these  "mistakes"  has  been  to  swing  the  pendulum  the  other  way  and  deny  fur¬ 
ther  permitting.  This  is  somewhat  the  case  for  residential  canal  develop¬ 
ment  in  the  State  of  Florida. 

This  process  was  examined  by  Snyder  (1976a,  1976b)  for  canal  developments 
in  Florida.  Early  engineering  criteria  were  (Figure  1 1 1 . 1 ) : 

1.  Navigable  depths  to  the  shoreline, 

2.  Maximized  front  footage  per  acre, 

3.  Increased  elevation  for  foundations, 

4.  Minimum  loss  of  property  to  water  area, 

5.  Rapid  drainage  of  rainfall,  and 

6.  Simplified  surveying  and  construction  methods. 

These  criteria  resulted  from  a  combination  of  trying  to  maximize  profits 
from  a  development ,  and  applying  hydraulic  practices  that  dated  back  a  num¬ 
ber  of  centuries.  In  many  cases,  the  criteria  were  borrowed  from  other 
types  of  development  and  caused  problems  in  tidal  canals. 


Figure  III.l.  Early  Canal  Design  (from  Snyder,  1976a) 


Some  of  the  results  of  the  above  approach  were: 

1.  Septic  tank  seepage, 

2.  Elimination  of  shallows, 

3.  Straight  canals  with  right  angle  bends, 

4.  Deep  narrow  canals, 

5.  Vertical  bulkheads, 

6.  Direct  stormwater  discharge,  and 

7.  Water  table  drawdown  from  local  wells. 

The  environmental  impact  of  canal  developments  using  the  above  criteria  are 
readily  visable  or  measurable.  They  include  domestic  waste  pollution,  the 
lack  of  natural  water  treatment,  poor  mixing,  limited  fish  and  animal  habi¬ 
tats,  stormwater  pollution,  and  saltwater  intrusion.  Because  of  the  result¬ 
ing  environmental  degradation,  most  regulatory  agencies  prohibit  the  con¬ 
struction  of  any  new  residential  canals  until  it  can  be  shown  that  such 
systems  are  compatible  with  the  site,  that  the  environment  is  not  degraded, 
that  all  regulatory  criteria  are  met,  and  that  the  system  is  rationally 
designed. 

The  regulatory  agencies'  response  to  the  problems  was  to  impose  criteria 
and  guidelines  that  they  believed  would  promote  more  environmentally  accep¬ 
table  designs.  Thus,  EP4  (1975)  recommended  that  "an  appropriate  canal 
depth  for  shallow  draft  pleasure  craft  should  be  no  more  than  four  to  six 
feet  below  mean  low  water,"  based  on  measurements  of  vertical  dissolved 
oxygen  profiles  and  numerical  flushing  models  (Figures  II. 7  and  II. 8). 
Other  criteria  specify  the  minimum  or  maximum  concentrations  of  given  water 
quality  parameters. 

Such  criteria  and  recommendations  were  a  response  to  the  problems,  and  some 
were  developed  through  research  and  observations.  Some  make  sense  in  envi¬ 
ronmental  terms,  and  others  need  interpretation  based  on  local  conditions 
and  other  factors,  As  an  example,  Horowitz  and  Bazel  (1977)  analyzed  the 
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process  by  which  an  advanced  wastewater  treatment  plant  was  proposed  for  a 
community  in  West  Florida.  They  found  that  the  justification  was  to  reduce 
phosphate  levels  below  regulatory  standards.  It  turned  out  that  the  ambient 
water  was  naturally  higher  in  phosphates  than  the  water  to  be  treated. 

This  is  not  an  attempt  to  criticize,  but  more  a  reminder  that  there  is  much 
to  learn  about  physical  processes  and  local  affects  before  a  rational  set 
of  design  criteria  can  be  imposed.  That  position  will  be  reached  by  trial 
and  error,  and  by  research.  For  example, the  EP4  Toxic  Standards  for  fresh 
water  are  now  being  revised  (Wu  et  al.,  198?)  because  in  some  areas  the 
animals  used  to  set  the  standard  do  not  exist.  In  canal  design  this  process 
of  learning  is  leading  to  a  new  approach  to  canal  development  (Snyder, 
1976a,  1976b;  Messenger  and  Reynolds,  1979).  Snyder,  for  example,  proposed 
a  new  set  of  design  criteria  (Figure  III. 2): 

1.  Offsite  sewage  treatment  or  package  plant, 

2.  Offsite  water  supply, 

3.  Retention  of  some  shallows, 

4.  Habitat  diversity--a  combination  of  vegetated  intertidal  shallows, 
sloping  rip-rap,  and  vertical  bulkheads, 

5.  Increase  circulation, 

6.  Promote  mixing,  and 

7.  Eliminate  scour  for  specified  conditions. 

These  are  not  a  rigid  set  of  criteria,  but  suggestions  that  should  lead  to 
more  environmentally  acceptable  designs  and  compliance  with  existing  regu¬ 
lations. Many  of  them  are  being  adopted,  and  many  are  being  further  investi¬ 
gated.  Van  de  Kreeke  and  Larsen  (1981)  and  Snyder  (1976),  for  example, 
discuss  water  quality  standards  for  marinas,  many  of  which  are  similar  to 
those  for  canals  and  have  arisen  because  of  the  same  concerns.  Howevpr,  in 
many  cases,  the  engineer  faces  the  task  of  selling  his  design  by  showing 
that  it  will  not  degrade  the  environment,  an  unenviable  task  when  there  is 
so  much  still  to  learn  about  what  constitutes  an  acceptable  design. 


II I -4 


"HF 


As  a  result,  the  design  engineer  must  select  an  evaluation  tool  that  is 
adequate  to  simulate  the  problem  being  investigated.  The  model  must  incor¬ 
porate  the  appropriate  physical  processes,  and  produce  results  that  a  per¬ 
mitting  agency  will  review  for  design  acceptability.  It  must  also  include 
features  that  are  elements  of  a  "rational"  design. 

If,  for  example,  wind  is  an  important  process  in  a  dead-end  canal  system, 
the  model  selected  should  be  capable  of  simulating  multilayer  flows  in  the 
vertical  direction,  as  this  process  cannot  be  simulated  using  the  area-mean 
velocity  of  a  one-dimensional  model. 

Permitting  agencies  often  want  to  know  "how  well  does  a  system  flush?"  or 
"what  will  be  the  minimum  DO  concentration  in  the  system?"  to  determine 
whether  a  proposed  system  will  conform  to  regulations.  They  may  wish  to 
know  what  the  sedimentation  rates  may  be,  and  whether  the  banks  will  be 
stable.  These  processes  should  be  included  in  the  model. 

Finally,  a  "good"  design  may  include  the  use  of  multiple  openings;  flow¬ 
through  systems;  large  upstream  water  bodies,  such  as  lakes,  to  flush  the 
system,  etc.  The  ability  to  include  such  features  must  be  built  into  the 
model . 

In  some  cases  the  combination  of  existing  processes  and  design  criteria  may 
simplify  the  modeling  approach.  For  example,  in  South  Florida  canals,  wind 
is  known  to  drive  circulation,  but  in  design  it  may  be  appropriate  to  ne¬ 
glect  its  effect  to  determine  whether  a  "worst  case,"  or  "poor  situation" 
meets  design  standards.*  In  such  a  case,  therefore,  it  would  be  appropriate 
to  use  a  numerical  model  without  a  wind  shear  term  to  make  this  evaluation. 

One  of  the  roles  of  numerical  models  is  the  establishment  of  rational  de¬ 
sign  criteria  through  model  sensitivity.  A  model  developed  by  Walton  (1978) 


*  Personal  Communication,  J.  van  de  Kreeke,  Rosenstiel  School  of  Marine 
and  Atmospheric  Science,  University  of  Miami,  Coral  Gables,  FL,  1982. 
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and  Walton  and  Christensen  (1980)  was  used  by  Morris  <1978),  Morris  et  al. 
(1978a,  1978b),  and  Hayes  (1982)  to  investigate  the  effect  of  canal  geome¬ 
tries  on  flushing  a  given  solute  distribution  or  inflow.  The  model  simu¬ 
lates  the  solute  transport  process,  and  allows  the  design  engineer  to  study 
transport  and  hydrodynamic  features  of  a  proposed  or  existing  system. 

Looking  at  tidal  forcing  only.  Figure  III. 3  shows  the  effects  of  increasing 
channel  length  and  inflow  rates.  Figure  1 1 1. 4  looks  at  the  effect  of  adding 
an  upstream  lake  to  improve  flushing  in  an  idealized  system  by  increasing 
the  tidal  basin  volume.  Such  studies  can  give  a  feel  for  system  response, 
but  Hayes  (1982)  tried  to  further  rationalize  the  process  by  developing  a 
series  of  design  curves  that  could  be  used  instead  of  the  computer  program. 

He  performed  a  dimensional  analysis  to  determine  a  series  of  dimensionless 
numbers  appropriate  for  canal  hydrodynamics  and  solute  transport  processes. 

By  performing  a  series  of  simulations  as  variations  of  a  standard  canal 
(Figure  III. 5),  he  developed  a  parametric  formulation  of  flushing: 

c  =  m(r,s)m1m2m3m4m5m6  (III  .1) 

where:  c  =  concentration  at  tidal  entrance, 

M(RS )  =  distribution  of  standard  canal,  and 
Mi -Mg  =  functions  of  dimensionless  numbers. 

W  1  odt.sxWsWe  lin'2’ 

time  to  maximum  concentration  at  tidal  entrance  (T), 
time  to  maximum  of  standard  canal  (T),  and 
functions  of  dimensionless  numbers. 

and 


where:  w  = 

D(R,S)  = 

VD6  = 
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CONCENTRATION  (ppm)  CONCENTRATION  (ppm) 


(L=LENGTH  OF  CANAL,  feet) 


DISTANCE  FROM  DEAD  END  (ft) 

Figure  III. 3.  One-Dimensional  Model  Sensitivity  Analysis 
(from  Morris  et  al.,  1978a) 
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Figure  III. 4.  Example  of  Model  Sensitivity  for  a  Canal  System 
Without  and  With  a  Lake  (from  Morris  et  al . ,  1978a) 


(ppm) 


M(R,  S) 


HYDRAULIC  RADIUS,  R  (feet) 

Figure  III. 5.  Standard  Runs  to  Determine  Concentrations  at  Tidal 
Entrance  of  a  Canal  System  (from  Hayes,  1982) 
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t10  .  N(R,S)NlN2N3N4N5N6 


(III. 3) 


where:  t^Q  =  time  to  flush  90%  of  solute  (T), 

N(R,S)  =  time  to  flush  90%  of  standard  canal  (T),  and 
N^-Ng  =  functions  of  dimensionless  numbers. 

Sensitivity  analyses  have  been  performed  on  marinas,  but,  to  date,  such 
studies  have  used  either  physical  models  or  field  measurements.  Nece  et  al  . 
(1975,  1976,  1979)  used  physical  models  to  study  geometric  effects  of  mar¬ 
ina  designs.  They  investigated: 

1.  Planform  geometry  aspect  ratio, 

2.  Ratio  of  entrance  cross-sectional  area  to  basin  planform  area, 

3.  Effects  of  rounding  corners  in  basin  interior, 

4.  Orientation  and  location  of  single  entrances,  and 

5.  Effect  of  two  entrances  versus  a  single  entrance. 

The  results  of  the  model  dye  studies  were  presented  in  terms  of  contours  of 
equal,  local,  per-cycle  exchange  coefficients,  E  (Figure  III. 6),  and  aver¬ 
age  tidal  flushing  (Figure  III. 7)  to  develop  some  recommendations  for  good 
marina  flushing. 

Slotta  and  Noble  (1977)  developed  a  statistical  model  based  on  observations 
of  sediments  in  existing  marinas  to  determine  flushing  potential  as  a  func¬ 
tion  of  geometric  properties.  From  the  study  they  developed  a  nomogram 
(Figure  III. 8)  for  use  in  the  design  of  marina  sitings. 

Such  studies,  using  numerical  models,  physical  models,  and  field  measure¬ 
ments  ,  are  a  sample  of  the  work  needed  to  develop  rational  design  criteria. 
In  the  future,  numerical  models,  as  they  perhaps  become  more  representative 
of  the  physical  and  water  quality  processes  existing  in  these  systems,  can 
be  used  to  simulate  many  design  alternatives  easily  and  quickly  to  learn 
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Figure  III. 7.  Sensitivity  of  Marina  Flushing  to 
Geometry  and  Depth  for  a  Single  Asymmetric 
Entrance  (Nece  et  al . ,  1979) 


Figure  I II. 8.  Nomogram  of  Acceptable  Marina  Goemetry  as  a 
Function  of  Planform  Area  (A),  Entrance  Width  (w),  and 
Entrance  Area  (a)  (from  Slotta  and  Noble,  1977) 
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how  such  systems  respond.  At  that  point,  they  will  be  useful  tools  to  de¬ 
sign,  evaluate,  and  establish  criteria  for  environmentally  acceptable  sys¬ 
tems  . 

To  conclude,  the  selection  of  a  model  for  evaluating  proposed  or  existing 
systems  must  include  the  consideration  of  both  physical  processes  and  de¬ 
sign  regulations.  The  engineer  must  be  aware  of  all  the  processes  affecting 
an  individual  system  before  a  decision  can  be  made  as  to  which  ones  to  in¬ 
clude  in  a  model.  This  decision  is  further  tempered  by  a  consideration  of 
what  issues  are  being  addressed  (such  as  worst  case  analyses),  what  pro¬ 
cesses  are  important  under  those  conditions,  and  what  may  be  neglected. 
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IV.  PAST  STUDIES  AND  EXISTING  MODELS 


Having  discussed  canal  and  marina  processes  and  design  considerations,  in 
this  chapter  we  will  review  some  past  studies  of  canal  and  marina  simula¬ 
tions,  and  examine  the  tools  used  for  analytic  evaluations.  We  will  begin 
by  considering  some  early  designs  before  looking  at  techniques  more  common¬ 
ly  used  today. 


IV. 1  Tidal  Prism  Analysis 


Most  of  the  early  designs  of  canal  and  marina  systems  were  based  on  a  sim¬ 
ple  tidal  flushing  analysis.  One  of  the  earliest  such  studies  was  by  Phelps 
and  Velz  (1933)  to  study  pollution  in  New  York  Harbor.  The  tidal  flushing 
analysis  is  based  on  volume  exchange  between  the  system  and  the  adjacent 
water  body.  Assuming  complete  mixing  at  each  exchange,  and  beginning  at  low 
tide,  the  flushing  can  be  expressed  in  terns  of  the  resulting  concentra¬ 
tion,  cN,  after  N  tidal  cycles: 


CN  "  co 


(  1  )N 
'T+F' 


where  cQ  =  initial  concentration  (dimensionless), 

N  =  number  of  tidal  cycles, 

F  =VW  f  3 

V  =  tidal  prism  volume  (L  ),  and 

=  volume  of  water  in  system  at  low  tide  (L  ). 


(IV. 1) 


or  by  considering  the  time  to  achieve  a  specified  dilution,  called  the 
flushing  time: 
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(IV. 2) 


t«  (c0/cN) 

*  ‘  Tn  d'  +  T)' 

where  N  =  number  of  cycles  to  desired  flushing. 

The  above  formulae  assume  a  one-dimensional  mean  exchange  between  the  sys¬ 
tem  and  adjacent  waters,  and  that  complete  mixing  occurs  on  each  tidal  cy¬ 
cle.  This  tends  to  overestimate  flushing  rates.  Modifications  were  made  to 
the  above  formulae  by  Gibson  (1959)  and  Ketchum  (1951)  who  divided  the  sys¬ 
tem  into  segments  and  assumed  exchange  between  each  segment  to  be  complete¬ 
ly  mixed.  A  second  modification  was  to  assume  some  incomplete  mixing  of  the 
tidal  exchange  volume.  This  rate  would  be  determined  empirically,  often 
about  0.8,  and  would  be  included  in  Equations  (IV. 1)  and  (IV. 2)  by  replac¬ 
ing  F  with  aF,  where  a  =  frational  rate  of  exchange. 

Even  today,  tidal  prism  analyses  are  used  to  evaluate  flushing  rates.  Kup- 
ferman  (1974)  developed  a  segmented  approximation  to  a  canal  system  and 
empirically  determined  exchange  coefficients.  In  studying  South  Beach  Ma¬ 
rina  in  Oregon,  Callaway  (1981)  used  a  simple  flushing  model  with  a  decay  of 
concentration  to  simulate  mixing  during  the  flood  tide.  His  results  (Figure 
IV.  1)  show  excellent  agreement  with  a  physical  model  of  the  system,  but 
also  show  that  the  method  can  tend  to  be  conservative  when  compared  with 
field  data.  Such  an  approach,  particularly  for  marinas,  does  appear  to  be  a 
reasonable  "back-of-the-envelope"  calculation  to  obtain  some  idea  of  the 
exchange  of  water  and  solute  mass  between  the  marina  and  the  adjacent 
waterway. 

It  is  further  interesting  to  note  that  such  an  approach  can  give  a  reason¬ 
able  estimate  of  numerical  model  results  of  flushing.  Lo  et  al .  (1976)  com¬ 
pared  a  simple  tidal  prism  method  with  a  one-dimensional  network  model  and 
achieved  quite  similar  results  (Figure  IV. 2). 


IV -2 


Relative  Dye  Concentration ,  Ci  /Co 


Hours  After  High  Water  Stack  Hours  After  High  Woter  Stock 


Figure  IV. 1.  Tidal  Prism  Flushing  Model  of  South  Beach  Marina 

( from  Callaway,  1981 ) 
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Figure  IV. 2.  Comparison  of  Tidal  Prism  and 
Numerical  Model  of  Flushing  (from 
Lo  et  al . ,  1976) 
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IV. 2  One-Dimensional  Models 


With  the  advent  of  computers,  the  design  engineer  was  able  to  use  more  so¬ 
phisticated  methods  to  evaluate  proposed  systems.  Instead  of  estimating 
flushing  rates  based  on  volume  exchanges,  hydrodynamic  and  mass  transport 
processes  could  be  directly  simulated  and  results  obtained  within  a  rela¬ 
tively  short  period  of  time.  The  models  allow  the  design  engineer  to  look 
at  not  only  flushing  rates,  but  also  to  evaluate  water  quality  conditions 
and  circulation  patterns. 

The  basic  governing  equations  in  one  dimension  are 
the  momentum  equation: 

iQu)  +  L  4  +  KB  iuiu  =  B*  +  *  (N,  (IV. 3) 

at  x  p  j*  11  Y  ax  l  ax 

the  continuity  equation: 


aA  aQ  _ 

iT+  7x 


*1 


and  the  mass  transport  equation: 


(IV. 4) 


where  Q 
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x 
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(Ac)  +  (Qc)  =  -V  (AE. 

3*  3  a  L 


3-f) 

a* 


-  KAc  +  qjCj 


flow  (L3/T), 

time  (T) , 

distance  (L), 

area  mean  velocity  (L/T), 

2 

acceleration  due  to  gravity  (L/T  ), 

2 

cross-sectional  area  (L  ), 
density  (M/L3), 


(IV. 5) 
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P  =  pressure  (M/LT2), 

k  =  friction  coefficient  (dimensionless), 

B  =  top  width  (L) , 

o  o 

4>  =  wind  shear  (L  /T  ), 

2 

N,  =  longitudinal  momentum  transfer  coefficient  (L  /T), 

L  2 
=  lateral  inflow  rate  per  unit  length  (L  /T), 

c  =  concentration  (dimensionless), 

=  longitudinal  dispersion  coefficient  (L  /T), 

K  =  decay  rate  (1/T),  and 

Cj  =  lateral  inflow  concentration  (dimensionless). 

Most  existing  one-dimensional  models  are  based  on  the  above  equations  or 
simplified  versions  of  them. 

Michel  (1973a,  1973b,  1973c)  used  a  model  with  simplified  versions  of  Equa¬ 
tions  (IV. 3)  and  (IV. 4)  to  study  canal  and  marina  circulation  in  Florida. 
He  assumed  that  the  velocity,  u,  could  be  given  by  Manning's  equation  in 
the  form: 


u  =  U486  r2/3s1/2  (IV.6) 

1  /3 

where  n  =  Manning's  roughness  coefficient  (T/L  ), 

R  =  hydraulic  radius  (L),  and 
S  =  hydraulic  gradient  (dimensionless). 

Equation  (IV. 6)  can  be  formed  from  Equation  (IV. 3)  by  assuming  a  balance 
between  the  pressure  force  and  the  frictional  resistance.  He  used  the  model 
to  look  at  the  changes  in  flows  resulting  from  proposed  geometric  changes 
due  to  marina  development. 

Another  simplification  of  canal  hydrodynamics  involves  the  assumption  of  a 
horizontal  water  surface.  As  was  noted  in  Section  II. I,  flows  in  canals  ex¬ 
hibit  very  small  water  surface  slopes  and  velocities.  Tide  gages  placed  at 
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either  end  of  a  system  would,  in  many  cases,  show  little  difference  in  ele¬ 
vation.  If  we  assume  the  water  surface  to  be  a  horizontal  plane,  tidal 
velocities  can  be  calculated  as  closed  form  functions  of  upstream  geometry: 


where  A$ 
h 


u  = 


ns  dh 
7T  M 


area  of  water  surface  (l  ),  and 
tidal  depth  (L). 


(IV. 7) 


The  above  model  (Walton,  1976a,  1976b;  Walton  and  Christensen,  1977)  was 
compared  with  a  model  of  Harleman  and  Lee  (1969),  based  on  the  full  equa¬ 
tions  (IV. 3)  and  (IV. 4).  The  results  never  differed  by  more  than  five  per¬ 
cent,  even  over  relatively  long  canals.  Such  variations  are  well  within 
error  limits  in  measurements. 


Knowing  the  velocity  allows  the  longitudinal  dispersion  coefficient  to  be 
calculated  as  a  closed  form  function,  based  on  the  formulations  of  Taylor 
(1953,  1954),  Elder  (1959),  Walton  (1978),  and  Walton  and  Christensen, 
(1980): 


El  =  KLhu  +  E0  (IV. 8) 

where  K.  =  dimensionless  longitudinal  dispersion  coefficient,  and 

L  2 
Eq  =  background  longitudinal  dispersion  coefficient  (L  /T). 

Equations  (IV. 7)  and  (IV. 8)  form  the  basis  of  the  hydrodynamic  and  solute 
transport  models  of  Langley  (1976),  Lee  (1977),  Walton  (1978),  and  Walton 
and  Chistensen  (1980). 

Westrich  (1978)  used  a  different  approach  for  studying  marina-type  basins 
adjacent  to  an  unsteady  main  stream  current.  By  considering  the  physics  of 
the  mass  transfer  process,  he  developed  a  model  that  simulated  mass  trans¬ 
port  in  both  the  main  channel  and  the  basin,  and  empirically  defined 
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exchange  coefficients  between  them.  The  model  performed  well  against  exper¬ 
iments,  but  more  work  needs  to  be  done  to  compare  the  results  against  field 
data  and  to  determine  the  range  and  sensitivity  of  the  exchange  parameters. 

In  most  cases,  however,  one-dimensional  models  are  developed  based  on  ap¬ 
proximations  to  the  full  governing  equations  to  simulate  hydrodynamics, 
flushing,  or  water  quality  transport.  Generally,  hydrodynamic  models  were 
developed  to  look  at  flows  through  large  canal  systems,  such  as  inter- 
oceanic  canals,  or  by  interbay  canals  (Keulegan,  1966;  Sagar,  1969;  Harle- 
man  and  Lee,  1969;  Boyd  et  al.,  1973;  Gardner  and  Pritchard,  1974;  Johnson, 
1974;  van  de  Kreeke  and  Dean,  1975;  Rives  and  Pritchard,  1978;  Najarian  et 
al . ,  1980).  Such  models  were  used  to  determine  heads  and  velocities,  and  to 
estimate  net  flows. 

When  evaluating  residential  canals  and  marinas  in  the  coastal  zone,  water 
quality,  rather  than  flow,  is  usually  the  issue.  In  fact,  as  relatively 
quiescent  systems  are  desired  for  navigation  and  boat  storage,  many  of 
these  systems  have  been  found  to  exhibit  poor  flushing,  and  have  water  qua¬ 
lity  problems  or  excessive  benthic  deposition. 

The  most  popular  one-dimensional  water  quality  model  is  based  on  the  link- 
node  concept  developed  by  Water  Resources  Engineers  in  the  mid-1960's  (Shu- 
binski  et  al.,  1965;  Feigner  and  Harris,  1970).  This  original  model  took 
two  different  directions,  toward  the  Dynamic  Estuary  Model  (DEM),  and  to¬ 
ward  the  RECEIV  model.  The  formulation  for  each  series  of  models,  however, 
is  very  similar,  and  between  them  they  have  modeled  more  estuaries,  canals, 
and  marinas  than  any  other  model. 

Recently,  Jettmar  et  al .  (1980)  used  a  modified  version  of  DEM  to  study  a 
canal  system  in  southeast  Florida.  Figure  IV. 3  shows  the  results  of  a  com¬ 
parison  between  computer  results  and  field  measurements.  Morris  et  al . 
(1980)  used  a  version  of  the  model  to  study  flushing  properties  in  another 
south  Florida  canal. 
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Other  water  quality  models  have  been  used  with  varying  degrees  of  success. 
Russo  and  McQuivey  (1975)  applied  the  simple  water  quality  model,  QUAL-1, 
to  a  south  Florida  canal.  Figure  IV. 4  shows  that  they  had  some  difficulty 
in  obtaining  of  good  fit  between  computed  results  and  observations. 
Brandsma  et  al .  (1973)  and  lo  et  al .  (1976)  applied  one-dimensional  branch¬ 
ing  models  to  investigate  mass  transport  and  flushing  in  marinas. 

Perhaps  the  most  comprehensive  investigation  of  water  quality  modeling  in 
canals,  however,  was  performed  by  EPA  in  the  early  1970 ' s  (EPA,  1973,  1975; 
Barnwell  and  Cavinder,  1975).  They  studied  a  number  of  canal  systems  from 
Florida  to  North  Carolina,  performing  dye,  water  quality,  and  modeling 
studies.  More  than  any  other  study,  this  one  allows  us  to  investigate  whe¬ 
ther  such  an  approach  is  appropriate  for  these  systems. 

Big  Pine  Key  Canal  is  one  of  several  canals  in  the  Florida  Keys  studied  by 
the  EPA  (1973,  1975).  The  canal  (Figure  IV. 5)  is  almost  1,600  feet  long,  40 
feet  wide,  and  10  feet  deep  at  mid- tide.  It  has  a  rectangular,  prismatic 
cross  section  and  is  aligned  with  its  longitudinal  x-axis  (defined  as  posi¬ 
tive  from  the  dead  end)  predominately  NNW. 

On  3  November  1973,  at  high  water  slack,  500  ml  of  Rhodamine  WT  dye  was  in¬ 
jected  at  mid-depth,  50  feet  from  the  dead  end  of  the  canal.  At  regular 
intervals  after  the  injection,  dye  concentrations  were  measured  using  a 
fluorometer  sampling  at  mid-depth  along  the  longitudinal  center  line  of  the 
canal . 

Once  the  data  were  collected,  EPA  tried  to  simulate  the  profiles,  using 
first  a  version  of  the  Storm  Water  Management  Model  (SWMM)  (Langer  et  al . , 
1971,  Metcalf  and  Eddy,  et  al.,  1971),  and  then  the  Columbia  River  Model 
(CRM)  (Callaway  et  al . ,  1969;  Callaway  and  Byram  1970).  The  reason  given 
for  turning  away  from  the  SWMM  model  was  that  (EPA,  1975): 
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Figure  IV. 5.  Longitudinal  Section  and  Cross  Section 
of  Big  Pine  Key  Canal  (from  EPA,  1973) 
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SWMM  did  not  satisfactorily  reproduce  observed  profiles.  The  field 
study  showed  the  dye  mixing  rapidly  throughout  the  length  of  the 
canal  and  then  slowly  "bleeding"  out.  The  model,  however,  shows 
the  dye  cloud  centroid  remaining  in  one  place  and  slowly  flatten¬ 
ing  out,  similar  to  classical  diffusion  without  advection.  As  the 
SWMM  neglects  diffusion  in  its  transport  equation,  this  apparent 
diffusion  is  probably  due  to  "numerical  dispersion",  or  “pseudo¬ 
dispersion",  associated  with  the  numerical  methods  employed  in 
thi s  model . 

Turning  to  CRM,  which  incorporates  a  dispersion  term,  several  of  the  canals 
examined  were  modeled  and  the  results  presented.  However,  the  measured  pro¬ 
files  were  reproduced  by  making  an  adjustment  to  the  dispersion  term,  ef¬ 
fectively  dividing  it  by  five. 

Lee  (1977)  also  tried  to  model  Big  Pine  Key  Canal  using  a  one-dimensional 
model  which  corrects  for  numerical  dispersion.  However,  in  obtaining  a  fit 
between  computed  and  observed  profiles  he  had  to  increase  the  dispersion 
term  by  two  orders  of  magnitude  (Figure  IV.6). 

From  these  results,  two  conclusions  may  be  drawn.  Firstly,  the  CRM  suffers 
from  severe  numerical  dispersion  when  used  to  model  low  energy  canals  such 
as  those  along  the  Gulf  and  East  Atlantic  coasts.  This  is  not  really  too 
surprising,  as  the  CRM  and  SWMM  are  interrelated  models,  both  having  been 
developed  from  the  Dynamic  Estuary  Model  (DEM).  The  inclusion  of  a  disper¬ 
sive  term  in  the  CRM  cannot  be  expected  to  produce  any  better  results,  as 
the  main  source  of  error  has  not  been  addressed. 

Secondly,  as  stated  in  the  EPA  report  (1975),  "the  dye  cloud  centroid  re¬ 
maining  in  one  place"  must  be  expected  in  a  model  which,  due  to  its  one¬ 
dimensional  form,  can  only  model  the  velocity  field  induced  by  the  astro¬ 
nomical  tide.  In  the  Big  Pine  Key  area  of  Florida,  the  tidal  range  is  on 
the  order  of  1  foot.  A  simple  volumetric  tidal  prism  analysis  indicates 
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that  maximum  velocities  of  no  greater  than  0.01  foot  per  second  can  be 
expected,  and  then  only  at  the  tidal  entrance.  At  the  dead  end  of  the 
canal,  the  tide-induced  velocities  will  be  much  lower,  and  certainly  not 
enough  to  move  the  dye  cloud  any  appreciable  distance. 

In  a  system  which  is  directly  open  to  the  ocean,  and  which  does  not  appear 
to  be  influenced  by  density  currents,  the  main  forcing  function  producing 
circulation  in  the  canal  must  be  the  wind.  This  effect  cannot  be  modeled  in 
one  dimension,  although  many  people  apparently  have  tried.  The  altering  of 
dispersion  coefficients  to  match  prototype  conditions  using  an  ill-defined 
model  can  only  give  results  good  for  that  particular  case--and  no  other. 
Any  predictive  analysis  based  on  such  a  model  calibrated  in  this  manner 
must  be  at  best  questionable. 

To  illustrate  this  point,  the  three-dimenrional  model  CANNET3D  (Walton 
1978,  Walton  and  Christensen,  1980)  was  run  for  this  case.  Three  layers 
were  used  in  the  vertical  and  an  estimated  wind  speed,  from  monthly  aver¬ 
ages  at  Key  West,  was  imposed.  As  can  be  seen  from  Figure  IV. 6,  a  closer 
approximation  to  the  shape  of  the  distribution  results,  using  reasonable 
values  of  transfer  coefficients. 

Part  of  the  problem  was  that,  in  many  instances  in  the  past,  one-dimensional 
water  quality  models,  which  were  really  intended  for  other  systems,  have 
been  applied  to  canals  and  marinas.  Specific  problems  were: 

1.  incomplete  physics  modeled, 

2.  incomplete  processes  modeled,  and 

3.  models  had  excessive  numerical  dispersion. 

As  we  have  seen  in  Chapter  II,  canal  systems  in  the  vertical  sense,  and 
marinas  in  the  horizontal  sense,  exhibit  two-dimensional  flows  on  many  oc¬ 
casions.  Attempts  to  study  these  systems  using  one-dimensional  models  have 
meant  that  major  processes  may  have  been  ignored,  and  that  their  effects 


had  to  be  lumped  into  a  group  of  uncertainties,  called  dispersion  and 
momentum  transfer  coefficients.  These  unknown  processes  are  then  modeled 
using  a  gradient  transfer  mechanism  that  may  not  even  approximate  the 
ignored  physics.* 

A  second  problem  is  that  many  of  the  numerical  models  are  not  set  up  to 
accurately  control  numerical  errors.  When  modeling  the  one-dimensional  ad- 
vection  equation: 


+  u  ^  =  0  (IV. 9) 

St  sx 

for  example,  using  an  explicit  upwind  difference  scheme,  we  obtain  a  numer¬ 
ic 

ical  dispersion  term  in  the  form  (Molenkamp,  1968)  En  — ^  where: 

sx 


E_  =  | uf AX  (1-  |u|it)  (IV. 10) 

n  —2—  -J5J— 

2 

where  En  =  numerical  dispersion  coefficient  (L  /T), 
aX  =  distance  interval  (L),  and 
At  =  time  interval  (T). 

Figure  IV. 7a  shows  the  effect  of  numerical  dispersion  for  this  scheme.  If 
the  center  of  mass  is  also  included  as  a  solution  parameter,  more  accuracy 
is  gained  as  shown  in  Figure  IV. 7b.  Including  the  second  moment,  or  width, 
allows  the  square  wave  to  be  advected  without  loss  of  shape.  This  latter 
scheme  is  the  basis  for  the  canal  model  of  Walton  (1978)  and  Walton  and 
Christensen  (1980).  Other  schemes,  such  as  the  central  difference  explic¬ 
it,  are  known  to  be  unstable.  These  schemes  are  discussed  at  length  in 
Roache  (1972). 


*  Personal  Communication,  J.  van  de  Kreeke,  Rosenstiel  School  of  Marine 
and  Atmospheric  Science,  University  of  Miami,  Coral  Gables,  FL,  1982. 
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In  modeling  rivers,  such  a  scheme  may  not  produce  too  severe  an  error,  or 

numerical  damping,  as  longitudinal  dispersion  coefficients  may  be  on  the 
2 

order  of  100  ft  /sec  or  more,  and  may  dominate  the  error.  In  canals  and 

marinas  where  ax  =  100  ft,  at  =  100  sec  and  u  =  0.5  ft/ sec  is  reasonable, 
2 

en  =  12.5  ft  /sec.  The  longitudinal  dispersion  coefficient,  EL,  observed  in 
these  systems  may  be  5  ft2/sec.  Thus,  a  model  simulation  of  such  a  system 
using  the  above  conditions  would  produce  severe  numerical  damping,  and  the 
results  would  be  questionable. 

The  problem  of  severe  numerical  dispersion  can  be  overcome  to  some  extent 
by  developing  a  scheme  to  minimize  the  error  produced.  The  schemes  of  Lee 
(1977),  Walton  (1978),  Walton  and  Christensen  (1980),  and  Fischer  (1970, 
1972)  specifically  addressed  this  issue  for  canal  systems.  Other  solution 
techniques,  such  as  those  of  Hall  and  Chapman  (1982),  Leonard  (1979), 
Fischer  (1978),  and  Glass  and  Rodi  (1982),  could  be  incorporated  into  canal 
and  marina  models  to  reduce  numerical  errors  in  the  mass  transport  solu¬ 
tion. 

We  have  shown  in  this  section  that  one-dimensional  models  cannot  always  be 
used  to  describe  and  simulate  the  processes  observed  in  canals  and  marinas. 
In  some  cases,  such  as  flow-through  systems  or  where  the  tide  dominates, 
they  may  be  adequate.  In  others,  such  as  dead  end  systems  or  systems  that 
exhibit  two-dimensional,  or  more,  flow,  they  may  be  less  than  adequate  to 
simulate  prototype  processes.  These  considerations  must  be  carefully  weigh¬ 
ed  before  selecting  a  modeling  approach. 


IV. 3  Multidimensional  Models 

Multidimensional  models  for  canal  simulations  are  rare,  and  when  they  do 
exist,  the  second  dimension  modeled  is  the  vertical.  This  is  reasonable  in 
most  cases  because  canal  systems  are  relatively  narrow,  and  physical  pro¬ 
cesses  such  as  wind  and  density  currents  tend  to  produce  vertical  variabil¬ 
ity  or  flow  reversals. 
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Such  models  have  usually  been  developed  to  overcome  the  problems  identified 
in  the  previous  section.  Monopolis  (1978)  realized  that  density-driven  cir¬ 
culation  could  not  be  simulated  well  using  a  one-dimensional  approach.  To 
simulate  the  effect  of  freshwater  inflows  to  saltwater  canals  in  south 
Florida,  he  developed  an  analytic  model,  incorporating  both  the  longitud¬ 
inal  and  vertical  directions,  to  predict  velocities  and  salinities.  He 
noted  that  the  results  (Figure  IV. 8)  shown  at  the  surface  were  similar  in 
form  to  observations,  but  more  work  needs  to  be  done  to  determine  the  util¬ 
ity  of  the  model  against  field  data. 

Walton  (1978),  and  Walton  and  Christensen  (1980)  developed  a  pseudo- three- 
dimensional  model  using  a  "flexing"  cell  structure.  Using  a  horizontal 
water  surface  assumption,  velocities  and  dispersion  coefficients  may  be 
specified  as  closed-form  functions  of  geometric  properties  using  the  con¬ 
tinuity  equation  with  head  given  at  the  tidal  entrance.  A  method- of- second 
moments  technique  was  used  to  minimize  numerical  dispersion,  and  the  model 
was  applied  to  several  canal  systems  in  the  State  of  Florida. 

The  model  generally  performed  well,  and  calibration  parameters  tended  to  be 
fairly  uniform  over  widely  varying  systems.  This  led  to  the  belief  that,  by 
including  the  correct  physical  processes,  even  in  simplified  form,  a  model 
could  be  developed  that  was  fairly  independent  of  calibration  conditions, 
and  this  might  be  useful  as  a  design  tool.  Comparisons  with  field  data 
showed  that  the  simplified  velocity  field  corresponded  satisfactorily  with 
velocities  in  a  wind-driven  system  (Figure  IV. 9)  and  in  a  density-driven 
system  (Figure  IV. 10).  Finally,  dye  studies  were  simulated,  and  again  the 
results  (Figure  IV.ll)  showed  good  agreement,  although  much  room  for  im¬ 
provement  remains. 

The  model  does,  however,  have  some  drawbacks  that  limit  its  usefulness.  It 
was  developed  to  model  upward  branching  systems  and  cannot  handle  multi¬ 
entrance  or  flow-through  systems.  Unfortunately,  designers  now  realize  that 
many  of  the  situations  the  model  was  meant  to  handle  (such  as  dead  end 
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Figure  IV. 9.  Comparison  Between  Model  and  Observed 
Wind-  and  Tide-Driven  Flow 
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Figure  IV.11.  Results  of  Dye  Simulation  in  57  Acres  Canal  System,  Florida 


channels)  should  not  be  built.  Another  deficiency  is  the  treatment  of  den¬ 
sity  currents,  which  was  empirically  based  on  observations.  The  model  was, 
however,  a  useful  beginning  and  could  serve  as  a  basis  for  further  develop¬ 
ment. 

As  in  canal  systems,  multidimensional  models  of  marinas  are  rare.  Most 
studies  have  used  either  one-dimensional  approaches,  or  psuedo-one- 
dimensional,  link-node  models.  Here  again,  most  models  are  simplified,  or 
neglect  some  processes  over  others. 

A  number  of  people  have  modeled  the  problem  of  shear-induced  circulation  in 
a  marina  basin  adjacent  to  a  large  water  body.  Abbott  (1977)  numerically 
simulated  the  combination  of  a  shear  flow  and  an  oscillatory  tidal  flow 
(Figure  IV.  12)  and  found  that,  to  calibrate  the  model,  he  had  to  increase 
the  momentum  transfer  coefficients  in  the  vicinity  of  the  basin  entrance. 

Askren  (1979)  solved  the  problem  of  a  steady- state  shear  flow,  and  modeled 
sediment  transport  and  deposition  within  a  rectangular  marina  basin  (Figure 
IV. 13).  No  mass  transfer  was  allowed  at  the  basin  entrance,  which  would 
have  required  the  inclusion  and  calibration  of  such  parameters. 

Chiang  and  Lee  (1982),  Falconer  (1980),  and  others  have  solved  the  much 
larger  (in  scale)  problem  of  circulation  within  a  harbor  (Figures  IV. 14  and 
IV. 15),  but  when  the  problem  is  reduced  to  the  scale  of  a  marina,  a  major 
drawback  occurs.  One  could  attempt  to  solve  the  governing  equation  impli¬ 
citly.  However,  for  Courant  numbers  above  five,  Vreugdenhil  and  Voogt 
(1975)  have  noted  that  phase  errors  can  begin  to  show  up  in  the  velocity 
field  depending  upon  the  number  of  grid  points/wavelength  employed.  For  im¬ 
plicit  calculations  below  a  Courant  number  of  five,  it  is  difficult  to 
justify  the  increased  computation  complexity  involved  in  solving  the  ap¬ 
proximate  equations.  Further,  for  an  explicit  solution,  a  Courant-governed 
hydrodynamic  solution  might  require  a  time  step,  at,  on  the  order  of  only 
several  seconds.  This  is  not  an  attractive  proposition  for  modeling  flush¬ 
ing  rates  and  water  quality  events  for  an  extended  period  of  time. 
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Figure  IV. 12.  Numerical  Simulation  of  Shear- 
Induced  Flow  in  a  Basin  (from  Abbott,  1977, 
courtesy  of  Danish  Hydaulic  Institute/ 
Oregon  State  University) 


Figure  IV. 13.  Numerical  Simulation  of 
Shear-Induced  Sediment  Transport  in  a 
Basin  (from  Askren,  1979) 
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The  models  described  in  this  section  serve  to  demonstrate  that  simple 
multidimensional  models  can  be  developed  to  simulate  non-unidirectional 
flows.  In  most  cases,  such  models  have  been  simplified,  both  in  terms  of 
dimensionality  and  process  assumptions,  through  judgment  of  observed  field 
data.  It  would  be  useful  to  develop  these  models  by  the  systematic  reduc¬ 
tion  of  a  more  complex  model  to  quantify  the  effects  of  assumptions  made. 


IV. 4  Other  Approaches 

Up  to  fairly  recently,  physical  modeling  of  marinas,  using  tidal  forcing, 
has  generally  been  chosen  over  numerical  modeling.  Even  considering  the 
limitations  of  using  distorted  models  for  determining  diffusion  and  mixing 
in  tidal  flows  (Fischer  and  Holley,  1971),  they  probably  better  represent 
the  physics  than  numerial  models  at  present.  Nece  et  al .  (1972,  1974,  1975) 
at  the  University  of  Washington  have  used  physical  models  successfully  for 
some  time.  It  is  also  interesting  to  note  that  they  and  Callaway  (1981) 
state  that  flushing  time  produced  from  these  models  is  in  fairly  good 
agreement  with  conventional  tidal  prism  analyses. 

However,  Slotta  (1982)  points  out  that  even  such  models  can  be  very  sensi¬ 
tive  to  induced  circulations.  As  an  example,  he  notes  that  the  flushing 
time  of  a  model  marina  was  reduced  from  ten  to  three  tidal  cycles  by  the 
induced  circulation  introduced  when  dye  was  added.  It  might,  therefore,  be 
useful  to  study  what  effects  wind-  and  density- induced  circulations  might 
have  on  the  flushing  times  of  marinas.  These  effects  may  perhaps  be  more 
easily  studied  by  using  a  nianerical  model. 

A  completely  different  approach  for  studying  canal  systems  was  proposed  by 
Snyder  (1980)  and  Snyder  et  al .  (1981).  Snyder  investigated  the  analogy 
between  the  components  of  electricity  and  oscillatory  flows  in  canal  sys¬ 
tems.  The  purpose  of  the  study  was  to  improve  mixing  and  circulation  pro¬ 
cesses  in  canal  systems  by  investigating  flow-through  systems  or  systems 
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that  mfght  resonate  with  an  adjacent  water  body  if  correctly  tuned.  When 
applied  to  some  existing  systems,  the  method  appears  to  work  well.  New  sys¬ 
tems  have  yet  to  be  built  to  test  the  theory  further.  This  approach  may 
make  an  interesting  tool  for  first-cut  analysis  of  a  proposed  system,  but 
it  appears  to  have  some  shortcomings  in  studying  effects  that  may  cause 
vertical  flow  variability. 
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V.  FUTURE  DIRECTION 


V.l  Scope 

In  the  previous  chapters,  the  present  state  of  canal  and  marina  numerical 
modeling  has  been  discussed.  The  main  conclusions  from  the  analysis,  to 
this  point,  are  that  we  still  know  relatively  little  about  such  systems, 
and  that  our  numerical  tools  are  not  always  adequate  to  make  constructive 
evaluations  of  proposed  designs. 

The  understanding  of  the  canal  and  marina  design  process  needs  to  be  im¬ 
proved  in  three  areas: 

1.  Processes, 

2.  Modeling,  and 

3.  Design  criteria. 

The  following  is  a  discussion  of  the  recommended  direction  that  future  work 
should  take. 


V. 2  Process  Identification 


Previous  studies  of  both  hydrodynamic  and  water  quality  processes  in  canals 
and  marinas  have  been  fairly  limited.  The  nature  of  the  data  collected  has 
been  "spotty,"  with  little  effort  to  monitor  a  system  over  a  long  period  of 
time.  This  effort  needs  to  be  consolidated,  and  a  highly  concentrated  mea¬ 
surement  program  needs  to  be  developed  for  three  or  four  canals,  and  three 
or  four  marina  systems,  that  cover  the  spectrum  of  expected  situations.  For 
residential,  tidal  canals,  the  following  areas  are  good  candidates: 
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1.  East  Coast  of  Florida, 

2.  West  Coast  of  Florida, 

3.  Gulf  Coast  of  Mississippi,  Louisiana, or  Texas,  and 

4.  Coast  of  California. 

For  marinas,  the  following  areas  should  be  considered: 

1.  Puget  Sound, 

2.  Oregon  or  Washington  Coast, 

3.  Coast  of  California, 

4.  Lower  Atlantic  or  Gulf  Coasts,  and 

5.  Upper  Atlantic  Coast. 

In  evaluating  suitable  choices  for  further  detailed  studies,  the  range  of 
sites  should  exhibit  the  various  processes  we  wish  to  quantify: 

1.  Tidal  exchange, 

2.  Wind-driven  flows,  both  in  flow-through  and  dead  end  systems, 

3.  Density-driven  circulation, 

4.  Geometric  effects,  such  as  multiple  entrances  and  upstream  storage 
areas , 

5.  Internal  energy  distributions  for  mixing  processes, 

6.  Sedimentation,  and 

7.  Water  quality  processes. 

Each  of  the  selected  sites  should  be  continuously  monitored  over  the  long 
term  (perhaps  one  year  or  more)  to  observe  the  following  parameters: 

1.  Tidal  elevations  at  several  points  in  t.he  system, 

2.  Wind  speed  and  direction, 

3.  Current  speed  and  direction, 

4.  Sedimentation  rates  and  di stributions , 

5.  Freshwater  inflow  quantity  and  quality. 
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6.  Water  quality,  as  salinity,  DO,  temperature,  etc.,  and 

7.  Atmospheric  parameters,  such  as  radiation  and  dew  point  temperature 
to  determine  effects  and  rates  of  water  quality  and  evaporation 
1 osses . 

Each  system  should  be  studied  in  detail  by  a  team  of  scientists  and  engi¬ 
neers,  and  their  results  presented  both  in  report  form  and  through  open 
discussions  at  a  conference.  The  aim  of  this  reporting  procedure  would  be 
to  identify  the  important  processes  of  canals  and  marinas,  and  to  propose 
methods  of  analytical  and  numerical  simulations. 


V.3  Numerical  Modeling  Approaches 

Several  levels  of  effort  are  recommended  for  numerical  approaches  to  simu¬ 
late  canal  and  marina  hydrodynamics  and  transport  processes.  The  first  is 
the  development  of  tools  that  will  simulate  all  the  observed  processes. 
This  may  consist  of  modeling  hydrodynamics,  sedimentation,  and  water  qua¬ 
lity  separately,  but  all  models  should  be  compatible.  The  second  is  the 
refinement  and  simplification  of  these  tools  for  practicing  design  engi¬ 
neers  . 

The  development  of  the  "scientific"  models  of  the  first  group  must  include: 

1.  All  processes  observed  from  field  measurements, 

2.  Highly  accurate  hydrodynamic  and  mass  transport  numerical  schemes, 
and 

3.  Selection  of  implicit  or  explicit  approaches. 

As  we  have  seen,  many  of  the  more  common  models  used  for  canal  and  marina 
simulations  neither  represent  the  important  physics  nor  accurately  model 
those  processes  included.  The  aim  of  developing  this  set  of  highly  accurate 


"scientific"  models  is  to  be  able  to  perform  sensitivity  analyses  to  deter¬ 
mine  the  relative  importance  of  terms  and  processes,  and  to  see  how  well 
these  models  simulate  prototype  conditions.  By  starting  with  all  the  pro¬ 
cesses  observed,  the  less  important  terms  and  processes  may  be  eliminated 
and  the  model's  ability  to  simulate  prototype  conditions  noted.  The  ulti¬ 
mate  aim  of  this  approach,  and  a  study  of  design  requirements,  would  be  the 
development  of  an  accurate,  yet  conceptually  simple,  suite  of  models  for 
canals  and  marinas  that  can  be  used  by  practicing  design  engineers. 

It  was  noted  in  the  previous  chapter  that  process  identification  alone 
would  not  necessarily  constitute  a  good  numerical  model.  The  accuracy  of 
the  numerical  approximation  can  sometimes  dominate  the  reliability  or  use¬ 
fulness  of  the  model  developed.  This  is  particularly  true  of  the  mass 
transport  equation  in  low  energy  systems  such  as  tidal  canal  networks,  in 
which  numerical  dispersion  can  dominate  that  being  modeled.  Lee  (1977), 
Walton  (1978),  and  Walton  and  Christensen  (1980)  tried  to  overcome  this 
problem  using  highly  accurate  numerical  schemes.  Unfortunately,  their  work 
was  limited,  and  further  studies  are  needed. 

There  are  many  estuarine  models  that  could  form  the  basis  for  the  model 
development  here.  A  thorough  review  was  given  by  Hinwood  and  Wallis  (1975a, 
1975b).  More  recently,  there  has  been  quite  a  lot  of  attention  paid  to  ac¬ 
curately  modeling  the  transport  equation,  particularly  in  two-dimensional 
estuarine  type  flow.  Such  work  (e.g.  Fischer,  1978;  Hall  and  Chapman,  1982; 
Holly  and  Preissman,  1977;  Leonard,  1979;  and  Glass  and  Rodi ,  1981)  could 
form  the  basis  of  mass  transport  model  development  for  canals  and  marinas. 

In  developing  numerical  models,  some  consideration  should  be  given  to  in¬ 
vestigating  explicit  and  implicit  approaches.  As  a  general  rule,  explicit 
approaches  are  simpler  to  conceive  and  program.  However,  for  canals  and 
marinas,  this  approach  may  lead  to  simulation  requiring  very  small  time 
steps  and  large  amounts  of  computer  time.  Traditionally,  this  has  been 
overcome  by  using  an  implicit  formulation.  The  problem  here  is  that  this 
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approach  tends  to  produce  significant  phase  errors  when  Courant  numbers 
larger  than  five  are  used  (Vreugdenhil  and  Voogt,  1975).  With  this  limita¬ 
tion,  there  tends  to  be  a  balance  between  the  Courant  condition  of  explicit 
techniques  and  the  extra  computations  involved  in  solving  the  implicit 
model .  Some  studies  for  reducing  the  phase  error  at  high  Courant  numbers 
may  be  worthwhile. 

Physical  models  have  been  used  for  many  years  to  evaluate  marina  designs 
(Nece  and  Richey,  1972;  Nece  and  Knoll,  1974;  Nece  et  al.,  1975,  1979),  but 
none  are  known  for  canals.  It  would  be  very  useful  in  selecting  sites  for 
comprehensive  measurement  programs  in  marinas  to  consider  at  least  one  for 
which  a  physical  model  already  exists.  The  aim  would  be  to  compare  physical 
and  hydraulic  models  against  field  data  to  assess  their  relative  merits, 
and  to  determine  whether  either  method  is  adequate  to  describe  the  pro¬ 
cesses  observed. 

In  both  the  canal  and  marina  studies,  it  would  further  be  useful  to  select 
sites  that  have  previously  been  numerically  or  analytically  modeled  during 
either  the  design  or  post-design  stages.  Such  studies  would  be  performed  to 
look  at  the  adequacy  of  present  design  methods,  using  the  more  sophisti¬ 
cated  techniques. 

The  design  engineers  require  a  set  of  models  that  enables  them  to  evaluate 
proposed  designs.  They  are  not  interested  in  the  detailed  adequacy  of  a 
particular  technique  or  model  to  simulate  their  system,  but  rather  an  inex¬ 
pensive,  accurate,  and  reliable  program  to  answer  their  questions.  Design 
models  need  not  contain  the  detailed  physics  of  "small  effect"  processes, 
but  rather  are  meant  to  look  at  the  effect  of  dominant  or  design  processes 
in  their  systems.  Reduction  in  dimensions  may  also  be  appropriate.  For 
example,  there  is  a  basic  need  to  determine  how  well  the  simple  tidal  ex¬ 
change  models  determine  flushing  rates  in  canals  and  marinas:  are  they  use¬ 
ful  design  tools? 
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In  developing  his  canal  model,  which  was  intended  for  use  by  design  engi¬ 
neers,  Walton  (1978)  and  Walton  and  Christensen  (1980)  assumed  a  horizontal 
water  surface  to  allow  the  velocity  field  and  dispersion  coefficients  to  be 
given  as  closed- form  functions  of  system  geometry.  He  also  assumed  a  super¬ 
imposition  of  velocities  representing  tidal-,  wind-,  and  density-driven, 
and  secondary  circulations.  This  approach  appeared  to  be  reasonable,  based 
on  a  limited  number  of  current  measurements.  It  would,  however,  be  desir¬ 
able  to  see  how  these  types  of  assumptions  hold  for  long-term  observations, 
and  to  use  the  observations  to  develop  closed-form  functions  for  the  var¬ 
ious  types  of  forcing. 

The  design  process  may  call  for  the  omission  of  one  or  more  processes  to 
evaluate  the  system  response  to  "worst"  conditions.  As  an  example,  van  de 
Kreeke*  considers  that  density-driven  circulation  dominates  in  the  long 
term  during  sunnier  months.  In  this  situation,  wind  effects  may  be  omitted, 
as  they  are  generally  light  and  variable  during  this  period.  Thus,  under 
some  design  conditions,  some  processes  may  be  omitted  and  models  simplified 
for  design  evaluation  purposes. 

Present  water  quality  models  are  generally  considered  over-complicated  for 
design  use  in  canals  and  marinas.  Large-scale  water  quality  models  are  ap¬ 
propriate  for  the  group  of  models  classified  as  "scientific."  The  design 
engineer,  however,  usually  would  like  something  less  sophisticated,  but 
adequate  to  address  his  evaluation  needs.  Also,  as  a  general  rule,  the  more 
complicated  a  model  is  and  the  more  calibration  coefficients  are  included, 
the  harder  it  is  to  obtain  a  satisfactory  fit  between  numerical  results  and 
observations.  More  often  than  not,  this  degree  of  sophistication  is  unwar¬ 
ranted  in  design  tools,  and  it  is  recommended  that  some  attention  be  given 
to  determining  which  water  quality  processes  are  important  for  design  eval¬ 
uation. 


*  Personal  Communication,  J.  van  de  Kreeke,  Rosenstiel  School  of  Marine 
and  Atmospheric  Science,  University  of  Miami,  Coral  Galbes,  FL,  1982. 
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The  simplest  form  of  water  quality  model  is  the  solute  transport  model,  in 
which  a  constituent  transport  is  included  in  a  hydrodynamic  model  for  the 
purpose  of  looking  at  flushing  rates.  As  an  example,  the  hydrodynamic  por¬ 
tion  of  DEM  ( DYNHYD )  is  being  extended  (Walton  et  al . ,  1982)  to  look  at 
salt  and  dye  transport  in  estuaries.  Such  simple  models  are  often  useful 
for  performing  a  first-cut  analysis  on  circulation  and  flushing  rates  in  a 
proposed  system. 

A  final  note  is  the  consideration  of  using  micro-computers  as  an  instrument 
for  design  models.  As  these  small  machines  become  faster  and  have  more 
storage  capacity,  they  may  be  appropiate  for  small  design  models.  Such  a 
development  would  be  very  useful  for  the  engineer  in  a  small  design  firm 
with  limited  computer  facilities. 

V.4  Establishing  Design  Criteria 

A  useful  product  of  numerical  model  development,  and  a  distinct  advantage 
over  either  the  prototype  or  physical  models,  is  the  ability  to  easily  per¬ 
form  sensitivity  analyses  to  establish  a  set  of  design  criteria.  These 
would  be  useful  for  both  the  design  engineer  and  the  permitting  agency. 
Ranges  of  these  studies  would  be  determined  from  analyses  of  present  and 
reconmended  measurements  in  canal  and  marina  systems.  Limits  of  acceptabil¬ 
ity,  in  terms  of  poor  water  quality,  flushing  rates,  and  sedimentation, 
must  be  known  a  priori  before  quantifying  the  limit  of  geometric  parameters 
to  comply  with  these  standards. 

At  this  time,  many  of  these  standards  already  exist.  For  example.  Federal 
regulations  define  what  constitutes  acceptable  water  quality  in  these  sys¬ 
tems.  Present  standards  should  be  reviewed  and  revised  in  view  of  new  com¬ 
prehensive  measurements,  and  should  take  into  consideration  local  factors. 
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Once  design  standards  have  been  established,  a  comprehensive  sensitivity 
analysis  needs  to  be  performed  to  identify  both  the  geometric  limits  of  a 
given  system,  and  also  which  design  elements  are  useful  in  achieving  the 
goal  of  environmental  compatibility.  Both  the  design  engineering  and  per¬ 
mitting  agency  need  to  know: 

1.  How  long,  wide,  and  deep  can  a  system  be? 

2.  Under  what  conditions  are  side  slopes  stable? 

3.  What  types  of  inputs  can  be  permitted? 

4.  Should  channels  be  aligned  with  the  wind? 

5.  Should  multiple  entrances  be  used? 

6.  Should  large  upland  water  bodies,  such  as  lakes,  be  connected? 

7.  What  is  the  effect  of  connecting  freshwater  and  saltwater  bodies? 

The  results  of  studies  such  as  these,  which  have  been  performed  on  some 
types  of  canal  systems  by  Morris  et  al .  (1978a,  1978b),  and  Morris  (1979), 
will  lead  to  a  better  understanding  of  what  elements  constitute  a  good, 
environmentally  acceptable  design.  Perhaps  this  work  can  even  be  extended 
in  the  manner  of  Hayes  (1982),  to  develop  a  handbook  of  design  curves  that 
would  allow  canal  design  to  be  peformed  without  having  to  resort  to  the 
computer. 

The  computer  is  a  useful  tool  for  the  scientist,  design  engineer,  and  per¬ 
mitting  agency  because  it  can  be  configured  to  run  a  large  number  of  test 
cases  with  relative  ease  and  economy.  We  must  be  careful,  however,  and  re¬ 
member  that  the  information  output  is  only  as  good  as  that  input.  Numerical 
models  can  be  a  tremendous  asset  if  they  are  well  formulated,  highly  accu¬ 
rate,  and  carefully  proved. 
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VI.  SUMMARY  AND  CONCLUSIONS 


In  this  report,  a  review  of  modeling  hydrodynamics  and  solute  transport 
processes  in  tidal  canals  and  marinas  is  presented.  By  considering  the 
three  areas  of  physical  processes,  modeling,  and  design  criteria,  we  see 
that  much  still  needs  to  be  learned  to  enable  the  engineer  to  design  such 
systems  to  meet  a  rational  set  of  design  criteria.  In  evaluating  such  pro¬ 
posed  designs,  the  permitting  agencies  are  in  the  same  relative  state  of 
uncertainty  as  to  what  constitutes  a  "good,  environmentally  acceptable" 
design.  Much  work  remains  to  be  done  to  develop  the  science. 

Future  direction  should  include  the  selection  of  several  varied  existing 
canal  and  marina  sites,  and  a  program  to  study  them  to  understand  the  fun¬ 
damental  physics  and  water  quality  processes  involved.  The  next  step  should 
be  to  use  this  information  to  develop  numerical  models  that  can  model  the 
important  processes  observed,  and  a  thorough  validation  of  their  formula¬ 
tion.  These  models  should  be  developed  as  engineering,  as  well  as  scienti¬ 
fic,  tools  to  allow  easy  use  by  practicing  design  engineers.  Finally,  a 
program  of  selective  model  and  measurement  sensitivity  should  be  performed 
to  expand  our  knowledge  of  rational  design  criteria.  These  should  be  cri¬ 
teria  that  make  sense  to  both  the  design  engineer  and  the  permitting  agen¬ 
cies,  and  that  lead  to  a  framework  within  which  both  can  work. 
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